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ON THE DISCREPANCY BETWEEN CHANDRA AND XMM TEMPERATURE PROFILES FOR A1835 
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ABSTRACT 

This short technical note addresses a large discrepancy between the temperature profiles for the galaxy cluster 
A1835 derived by Schmidt et al. (2001) using Chandra and by Majerowicz et al. (2002) using XMM. The causes 
of this discrepancy may be instructive for the Chandra and XMM cluster analyses in general. The observation 
used by Schmidt et al. was affected by a mild background flare that could not be identified by the usual technique. 
This flare biased upwards the measured temperatures at large radii. The remaining discrepancy appears to be 
due to the XMM PSF scattering that was not taken into account in the published analyses. While the XMM PSF 
is narrow, the surface brightness of a typical cluster also declines very steeply with radius. For the moderately 
distant, cooling flow cluster A1835, about 1/3 of the observed XMM brightness at any radius is due to the PSF 
scattering from the smaller radii. As a result, the contamination from the bright cool cluster center biases low the 
measured temperatures near the core, and in general, any temperature gradients are underestimated. 
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1. INTRODUCTION 

A1835 is a moderately distant (z = 0.25), symmetric cluster 
with a strong central X-ray brightness peak. Using Chandra 
data, Schmidt et al. (2001) measured its temperature profile out 
to r = 4'. It revealed a cool center and a temperature increase 
to ~ 13 keV at r = 0.5- 1 h^ Mpc. However, Majerowicz et 
al. (2002) derived a very different temperature profile from the 
XMM data. While they confirm the presence of the cool central 
region, their outer temperature reaches only ~ 7 keV. The two 
profiles in the common range of radii are shown in Fig. [j]. Such 
a discrepancy at the qualitative level merits detailed investiga- 
tion. I first re-analyze the Chandra data. 



20 



' 5 



' 



Chandra, Schmidt et al. 



XMM, Majerowicz et al. 



1 Mpc 



r, arcmin 



FIG. 1 . — Published Chandra (Schmidt et al.) and XMM (Majerowicz et 
al.) temperature profiles for A1835. Errors are 90%. 



2. CHANDRA ANALYSIS 

There were 2 Chandra observations of A1835, with the clus- 
ter in the ACIS S3 chip in both. One (OBSID 495) with a 19 ks 
exposure was performed on 1999-12-11, and another (OBSID 
496) with an 1 1 ks exposure was done on 2000-04-29. Schmidt 
et al. used only OBSID 495 for their spectral analysis. They 
used 0.5-7 keV energy band for spectral fitting and the public 
blank field dataset for background modeling. In OBSID 495, 
only one backside-illuminated (BI) chip, S3, was used, and it 
was almost completely covered by the cluster emission. This 
presents a difficulty for identifying the background flares (see 
Markevitch 2001a) using the background light curve. A light 
curve extracted from a region of the S3 chip far from the clus- 
ter center does not vary beyond the normally acceptable ±20% 
range from the mean (in agreement with the Schmidt et al. con- 
clusions), although due to the cluster contamination, one can- 
not compare that mean value to the nominal background rate. 
In OBSID 496, one can also use S 1 (another BI chip) for flare 
detection; neither SI nor S3 (its part far from the cluster center) 
show any flares. 

Since the quiescent background is known to vary on long and 
short timescales at a few percent level, for each observation, I 
calculated the normalizations of the corresponding blank-field 
datasets from the ratios of count rates in the 10-12 keV band 
(which is free of the source emission and of the most common 
variety of flares). As expected, such normalizations were within 
a few % of the respective exposure ratios. 

ARFs and RMFs are calculated as described in, e.g., Marke- 
vitch & Vikhlinin (2001). The recently discovered low-energy 
ACIS quantum efficiency change was corrected; it is not impor- 
tant for the present analysis. Nh was fixed at the Galactic value 
and abundances at 0.3 solar. 

For OBSID 495, using the same background and energy band 
as Schmidt et al., I obtain a temperature profile very similar to 
theirs. 



2.1. Soft Galactic background excess 
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FIG. 2. — Chandra A1835 temperature profiles from observations OBSID 
495 (the one analyzed by Schmidt et al.) and OBSID 496, derived using the 
0.8-8 keV band. The soft Galactic background excess component was ac- 
counted for. Errors are 90%. 
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FIG. 3. — Chandra radial brightness profiles from OBSIDs 495 and 496 in 
the 2.5-7 keV band (most sensitive to the background flares). The Sx scale 
is arbitrary but the observations are normalized by their exposures. Errors are 
icr. There is a clear background excess in OBSID 495. 



However, as pointed out by Majerowicz et al., A1835 is pro- 
jected onto the Galactic North Polar Spur (Snowden et al. 1997) 
and thus has anomalously high sky background at E <^ 1 keV. 
Since the public ACIS background datasets represent the aver- 
age high Galactic latitude background (see Markevitch 2001b), 
this anomaly should be taken into account. To evaluate this 
soft background excess, I extracted spectra from ACIS chips 
far from the cluster, subtracted the nominal background, and 
fit the residuals by an arbitrary model (I used low-temperature 
MEKAL with solar abundances). This best-fit model (renor- 
malized by the ratio of the solid angles) can then be added to 
the fits in the interesting cluster regions. This properly takes 
vignetting into account, assuming that the excess component is 
indeed celestial in origin and does not vary on the ACIS FOV 
scale. This method was used for Chandra analysis by Marke- 
vitch & Vikhlinin (2001), and a similar method (without the 
modeling step) was used for XMM by Pratt et al. (2001) and 
Majerowicz et al. 

For this procedure, a region of chips 12 and 13 at r > 10' 
(3 Mpc) was used for OBSID 495, and the SI chip was used 
for OBSID 496. I created subsets of the blank-sky background 
datasets for the respective time periods ('B' for OBSID 495 and 
'C for OBSID 496) that included exactly the same blank fields 
for S3 and those other chips (the next release of the public back- 
ground files will include such subsets). For the 'C period, such 
a subset unfortunately is short and dominated by two pointings 
toward the Galactic Spur. This resulted in negative soft excess 
for OBSID 496, which was fit by a negative model. In OBSID 
495, there is a real excess above the blank field background. 

The best-fit soft excess (or deficit) was included in the cluster 
fits, and a low energy cutoff of 0.8 keV was adopted to mini- 
mize the associated uncertainty. The resulting radial profiles for 
the two observations are shown in Fig. ||. The two are inconsis- 
tent, and OBSID 495 clearly has a problem. 

2.2. Flare background correction 













-h 









channel energy (keV) 



FIG. 4. — Spectrum from the region of the ACIS S3 chip outside r = 4.5' 
(1.4/jjg Mpc) from the cluster center, with the nominal background subtracted. 
There is a clear unphysical hard excess, which in the 2.5-7 keV band is con- 
sistent with the background flare model (shown by histogram; see text). The 
residuals at lower energies are due to cluster emission and the Galactic soft 
background excess. 

The most obvious candidate for the problem in OBSID 495 
is a mild background flare, which are often observed in the BI 
chips. It may be too faint or too constant to be identified in the 
light curve of the relatively short observation (especially if one 
is using a region containing cluster emission in addition to the 
background). Spectral shapes of the flare and quiescent back- 
ground components are such that the 2.5-7 keV band is most 
sensitive to flares. Fig. || shows brightness profiles in this en- 
ergy band for the two observations, after subtracting the nomi- 
nal background. OBSID 495 shows a clear background excess. 

Figure [| shows a spectrum of this excess, extracted from the 
region of the S3 chip outside r = 270" from the cluster center. If 
the cluster brightness profile from OBSID 496 in Fig. || is cor- 
rect, the true cluster brightness in this band at these large radii 
is a small fraction of the background excess. 



A1835 CHANDRA-XMM DISCREPANCY 
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FIG. 5. — Chandra profiles from OBSIDs 495 and 496, now with the cor- 
rected background flare component in OBSID 495. Errors are 90%. 



10 



Chandra corrected 



XMM 



r, arcmin 



1 hr' Mpc 



FIG. 6. — The corrected Chandra temperature profile (average of 2 observa- 
tions) and the XMM profile. Errors are 90%. 



There are at least two species of the background flares ob- 
served in the BI chips. The most frequent one affects only the 
BI chips and is not seen in FI chips. In OBSID 495, there are 
no flares in the FI chips (this can be established confidently us- 
ing the regions free of cluster emission covered by those chips). 
Thus the putative flare in chip S3 should belong to this species. 
Such flares appear to have the same spectral shape regardless of 
their brightness, and are well-described by a power law with the 
photon index 7 w -0.1 and an exponential cutoff at E « 5 keV 
(Markevitch et al. 2002a) without the application of the tele- 
scope or CCD efficiency ('arf none' in XSPEC). Figure |] shows 
this model with the normalization fitted to the spectrum in the 
2.5-7 keV band. The excess at E > 2.5 keV is fully consistent 
with this model; the residuals at lower energies should be due 
to the soft Galactic excess (uncorrected for this exercise) and 
the cluster emission. At E = 3 - 5 ke V, the flare brightness is 
30-40% of that of the nominal background. 

Thus, the hard excess in OBSID 495 is consistent with a mild 
residual background flare. We can assume (somewhat arbitrar- 
ily) that this component is spatially uniform and subtract it from 
the cluster spectra, normalizing it according to the solid angles. 
This method was used in Markevitch et al. (2002b). The result- 
ing corrected temperature profile for OBSID 495 is shown in 
Fig. |] along with the original profile for the unaffected OBSID 
496. A ±40% uncertainty (90% confidence) for the flare com- 
ponent normalization is included in quadrature. The effect of 
the flare correction is significant at r > 1'. After the correction, 
the two Chandra observations become more or less consistent. 
(Note that the original Schmidt et al. profile shown in Fig. 1 is 
consistent with the corrected profile; the competing effects of 
the soft Galactic background excess and the hard flare excess 
partially compensated for each other in that work.) 

Figure || overlays the corrected Chandra temperature profile 
on the XMM profile; there is still a significant discrepancy. 
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3. XMM PSF EFFECT 

'http://heasarc.gsfc.nasa.gov/docs/xmm/uhb/xmm_uhb.html, sections "X-ray telescopes", "On-axis PSF' 
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FIG. 7. — Fraction of the observed flux in each XMM annulus originating 
from the same annulus in the sky, and fractions scattered by the telescope 
from the central circle and from all inner regions including the center. 

The obvious candidate for the remaining discrepancy is the 
XMM PSF, which was not included in any currently published 
XMM analyses. For a quick estimate of the PSF effect, I use 
on-axis MOS1 PSF in the 0.75-2.25 keV band presented in 
"XMM-Newton Users' Handbook". 1 It can be described by the 
function Sx oc (l + r 2 /r^)~ a with r c = 3.5", a = 1.36. Ignor- 
ing its dependence on coordinate and energy, one can estimate 
what fraction of the observed XMM flux in each region origi- 
nates in the same region in the sky and how much is scattered 
from other regions. I used the Chandra image of A1835 as a 
true brightness model (the Chandra PSF can be ignored for the 
present purpose). The model image was cut into annuli (same 
as used by Majerowicz et al.), each annulus is convolved with 
the PSF, and contributions from each annulus in the sky into 
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each annulus in the image were calculated. 

The result for the inner 5 annuli (those completely covered by 
Chandra) is shown in Fig. ^[ Even though the PSF 90% encir- 
cled energy radius is only 45-50", the cluster brightness profile 
declines so steeply with radius that in each annulus, there is a 
~ 10% or higher contribution from the central peak. Further- 
more, every annulus has a ~ 30% contamination from the inner 
regions, with the obvious effect of smearing any temperature 
gradients in the cluster. 

A simple XSPEC simulation (using EPIC responses and the 
exposure time, energy band and spectral binning given by Ma- 
jerowicz et al.) shows that if the temperatures from the Chandra 
profile are mixed in the proportion shown in Fig. one can ob- 
tain single-temperature EPIC fits very close to those reported. 
In particular, for the 2nd annulus, I obtain 6.4-6.8 keV com- 
pared to Chandra's 8 keV. For the 3rd annulus, I obtain 7.0-7.7 
keV compared to Chandra's 9 keV. 

4. LESSONS 



1. When ACIS background is critical, one should ensure 
that the quiescent background rate is consistent with the nom- 
inal rate (the one in the background datasets), even if there 
are no obvious signs of background flares in the background 
light curve. For ACIS BI chips, the 2.5-7 keV band is best for 
this purpose because it is sensitive to the most common flare 
species. The next version of the public background datasets for 
the BI chips will be filtered using this energy band (instead of 
the presently used wide band). 

2. The XMM PSF has to be taken into account for the analy- 
sis of clusters with peaked X-ray brightness profiles. (Monique 
Arnaud communicates that it will be included in their forthcom- 
ing reanalysis of the A1835 data.) 

Incorrect modeling of the Chandra background and disre- 
gard of the XMM PSF appear to explain most if not all of the 
discrepancy between the recently published A1835 temperature 
profiles. 

This analysis was supported by NASA contract NAS8-39073 
and grant NAG5-9945. 



REFERENCES 



Mayrnwir7, S , Npupiann, D. M., & Reiprich, T. H. 2002, A&A, submitted 



M. 



M. 



urn: 

istro-ph/0202347| ) 

E53EE M I a mpmn "APTS hnrlrgrniinH" 



http://cxc.harvard.edu/contrib/maxim/bg 



http://cxc.harvard.edu/contrib/maxim/acisbg/data/README 



Markevitch, M., & Vikhlinin, A. 2001, ApJ, 563, 95 

Markevitch, M., et al. 2002a, in preparation (^'Chandra spectra of the soft X-ray 

diffuse background") 
Markevitch, M., et al. 2002b, ApJ, 567, L27 

Schmidt, R. W., Allen, S. W., & Fabian, A. C. 2001, MNRAS, 327, 1057 
Snowden, S. L., et al. 1997, ApJ, 485, 125 



